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ABSTRACT 

We propose the B c — > B* / y decay as the most suitable probe for the flavour changing neu- 
tral transition c — > wy. We estimate the short and long distance contributions to this decay 
within the standard model and we find them to be comparable; this is in contrast to radia- 
tive decays of D mesons, that are completely dominated by the long distance contributions. 
Since the c — > transition is very sensitive to the physics beyond the standard model, the 
standard model prediction Br(B c — > B**y) ~ 10~ 8 obtained here opens a new window for 
future experiments. The detection of B c — ► B*^ decay at branching ratio well above 10 -8 
would signal new physics. In addition we study the s — ► d^y transition in B s — > B^ decay 
and we find it to be dominated by the long distance contribution. We use the Isgur-Scora- 
Grinstein-Wise (ISGW) constituent quark model for the calculation of these decays. 



1. INTRODUCTION 

Flavour changing neutral current (FCNC) transitions occur in the standard model only 
at the loop level. Hence, they are very rare in the standard model and they present a suitable 
probe for new physics. The FCNC transitions in the down-quark sector are relatively frequent 
due to the large mass of the top quark running in the loop and the transition b — > s has 
indeed been observed Jl| . The FCNC transitions in the up-quark sector are especially rare in 
the standard model due to the small masses of the intermediate down-like quarks that run 
in the loop. For these transitions, the standard model represents a small background for the 
possible contributions arising from some new physics. At present, only upper experimental 
limits on the FCNC transitions in the up-quark sector are available [0]. 
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In the present paper we study the transition c — > wy, which is the most probable FCNC 
transition in the up-quark sector within the standard model. To probe the c — > wy transi- 
tion we propose the radiative beauty-conserving decay B c — > B^y; the B c meson has been 
detected recently at Fermilab ||. We estimate the short distance (SD) and long distance 
(LD) contributions to B c — > B*j decay within the standard model and find them to be 
comparable, which allows us in principle to probe c —>■ wy transition in this decay This 
is in contrast to previously discussed D meson decays, where the dynamics is completely 
dominated by the LD contributions [4-7] and it is impossible to extract the short distance 
c — > wy contribution from the experiment. 

In addition we study s — > d'y transition in B s — > B^y decay. The dynamics of B s — > B^y 
is very similar to the dynamics of B c — > B*^ due to the spectator b in both decays. In 
contrast to B c — ► B*j, the decay B s — > B^y is found to be dominated by LD contributions 
and consequently the signal of new physics is not expected in this decay. 

In Sec. 2 the B c — > B*y decay is studied: we define SD and LD contributions, present 
the model and the results. The same reasoning is applied to B s — > B* d y decay is Sec. 3. We 
conclude with a summary in Sec. 4. 

2. B c -f B* u y DECAY AND c -> m 7 TRANSITION 

2a. The short distance contribution 

The SD contribution in B c — > B*j decay is driven by FCNC c — > wy transition and b is a 
spectator. The c — > wy transition is strongly GIM suppressed at one-loop, QCD logarithms 
enhance the amplitude by two orders of magnitude ||, while the complete 2- loop QCD 
corrections further increase the amplitude by two orders of magnitude 0. The Lagrangian 
that induces the c — > wy transition is given by 

C-Td" = -^^ V -Ks <r m (») ua^\m}-±^ + mJ-^-]c F, v . (1) 

The appropriate scale for C7~ > " 7 (/i) in B c — > i?*7 decay is /x = m c (and not \i = nib), since b 
is merely a spectator in the SD process. The 2-loop QCD calculation was performed in [|J, 
giving C 7^" 7 (m c ) = -0.0068 - 0.020i 

The corresponding amplitude for B c — > B*j(q, e) decay is proportional to e*^q u {B*\ua^ u (1± 
7 5 )c|5 c ) taken at q 2 = 0, which can be expressed in terms of the form factors i*i(0) and F 2 (0) 

e;(B* u {p',e')\uia^q u c\B c {p)) g , =0 = ie^e^p^O) , 
e^^^eOlm^g^cl^W)^ = [(m% c - m 2 B «)e* • e*' - 2(e*' • g)(p • e*)]F 2 (0). (2) 



The form factors defined in Eqs. (@) will be calculated using the ISGW model |T(J later on 
in this section. 

2b. The long distance contributions 

The long distance contributions will be calculated by using the nonleptonic weak La- 
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grangian, which on the quark level can be written as || 

G 

C eff (Ac = 1) = -^V t ^.V^.[a 1 (uq i )^_ A (q j c) v ^ A>lx + (hfa)v-A(Wi)v-A,p] , (3) 

where (ipi?p2)v-A = ^i7 M (l — 75)^2, Qi, Qj are the down quarks d, s, b and ai, 02 include the 
QCD corrections . 

Quite generally, the LD contributions to B c — > B*'-/ decay can be separated into two 
classes related to the two terms of (0), as performed previously H for D —> decays. The 
class (I) , called also the VMD contribution, is related to the a 2 term (H) and corresponds to 
the processes c — > uqiqi followed by — > 7, while b is the spectator in B c — > B*j decay. 
At the hadron level the q\qi ~^ 7 transition is expressed using the vector meson dominance 
(VMD) and the corresponding diagram is depicted in Fig. (la). The class (II), called also the 
pole contribution, is related to the ai term (H) and corresponds to the process cb — > ub with 
the photon attached to incoming or outgoing quark lines. Selecting the lowest contributing 
states, the pole contributions are depicted in Fig. (lb). 

We turn now to the estimation of these two classes of contributions and we start with 
the VMD contribution (class (I)) represented by Fig. (la). The underlying quark processes 
are c — > uss(dd) with ss, dd hadronizing into vector mesons <j), p, to which then turn to a 
photon, while b remains a spectator. We neglect the contribution of bb —>■ 7 in view of the 
large mass of T. The relevant part of the Lagrangian, after using the relations among CKM 
matrix elements, is 

C$ = ~^aMV cs V: s uY(l ~ 7s)c [*y M (l - 7s)* - d 7 „(l - l,)d] . (4) 

The appropriate scale for 02 (/i) in B c — > 5*7 decay is n — m c , since b is again merely 
a spectator in VMD contribution. Thus, we may use a2(m c ) = —0.5, as obtained in the 
successful phenomenological fit to D meson decays fTl"] . Defining (V(q, e)|V^|0) = gy(g 2 )e 



and using the factorization approximation, the effective Lagrangian that induces the VMD 
contribution is given by 

£ C VMD e) = -^j=a 2 (m c )V cs V* s C' VMD u^{l - 7 5 )c e* , (5) 

were 

is obtained by assuming gy(my) = <3v(0), with the mean value and the error in @ calculated 
from the experimental data on T{V — > e + e~) 0. Note here the remarkable GIM cancellation 
carried over to the hadronic level. 

Lagrangian fl5|) implies that the VMD amplitude for B c — > 5*7 (g, e) is proportional to 
e* (5* 1^7^(1 — 75)c|i? c ) taken at g 2 = 0. For the hadronic matrix elements, one defines 
appropriate form factors for the vector and axial transitions as follows f9|: 

2i 

(B:(p>, e') \u^(l - J 5 )c\B c (p)) = e^ey p 7 V(q 2 ) + (m Bc + m B ,)e^A 1 (q 2 ) 

TTT'Bc + m B* 

e^_q_ {p + p y A 2(g2) _ 2 f^^^ _ Ao(g2)] (?) 

ws c + m B * g 2 
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The requirements of the finite matrix elements at q 2 = [ I I ] and of gauge invariance lead 
to the relations among the various form factors ||, which imply A (0) = A 3 (0) = and 
^2(0) = [ijn Bc + m B *) I (m Bc — m B *)]Ai(0). The same relations are obtained by using the 
prescription that the photon couples only to the transverse polarization of the current || |P2| . 
Accordingly, the VMD amplitude will be expressed in terms of two form factors only, V(0) 
and Ax(0). 

At this point, we remark that the form factors Fi, F 2 , V and Ai, needed for the SD 
and VMD amplitudes cannot be safely related using the Isgur-Wise relations fL3 |, since the 



masses of b and c quarks composing B c meson do not permit the b quark to be at rest. 
Therefore we shall determine the corresponding form factors at q 2 = independently, using 



the ISGW model 10 



We now turn to the discussion of the LD contributions of class (II) , the pole contribution, 
where the quark process cb — > ub is driven by 



and the photon line is attached to any of four quark lines. In terms of hadronic degrees 
of freedom this diagram is given in Fig. (lb), where the white box represents the action 
of the Lagrangian (§) (we have neglected the contribution of the scalar and axial poles). 
Considering the scale for ai(/i) in cb — > ub, it is difficult to decide between /1 = m c or fi = m b , 
since b is not spectator in the pole contribution. As the difference between a\{m c ) = 1.2 and 



ai(rrib) = 1.1 [|ll[ an d is not essential, we take ai(rrib) = 1.1. Note that the pole contribution 
is relatively small due to the factor V c bV* b in (^). In D meson decays, the corresponding 
factor V CS V* S is much bigger, which makes the pole contribution dominant over the SD and 
VMD ones [f|, [5], |6|. The different CKM matrix elements in the pole contributions of B c 
and D decays is a major factor in establishing the B c — > B*"f decay as more suitable for the 
investigation of c — > wy than the D decays. 

To evaluate the amplitude for the pole diagrams given in Fig. (lb) we define 

(0\A^\P) = f P p, (9) 
(V\V,\0) = g v e; 
A(P(p) - V (p', e')j(e)) = Upee^eyMp , 

where fi Bc , fi Bu , f Bc , f Bu , g B * and g B * will be determined using ISGW model. 

2c. The amplitude 

Using the above Lagrangians and form factor decomposition of Eqs. (0), (|7|), (|9]), the 
final amplitude for B c — > B*'-/ containing SD and LD contributions can be expressed as 

A(B c (p) - Blip', e'h(q, e)) = ie*el[A PV (p»p» - g^p ■ q) + iA PC e^p' a p p ] , (10) 

where 



Apv = -^e[V cs V: d 



" {i ' K) <m c -m u )F 2 {U) + Za 2 {m c )C' VMD - Al{0) 



n 9 v ^ "/ *\ / 1 £ \ v 1V1 U 

2ir 2 m Bc - m Bt 
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PC 



v/2 

+ V cb V* b ai{m b 



4tt 2 
^b c 9b* c 9b* 



(m c + m u )Fi(0) + 2a 2 (m c )C[ 



V(0) 



VMD , 



m Bc + m B * 



+ 



VB u m B JbJb u 



m n — m 



B u 



(11) 



The first term in Eqs. (|Tl|) comes from SD contribution, the second term from VMD con- 
tribution and the third term from the pole contribution. The decay width is then given 
by 



1 / rn? Bc - m\ 



1-4 



PV 



An \ 2m Bc 

2d. The model 



+ \A 



PC\ 



(12) 



To account for the nonperturbative dynamics within the mesons we use the nonrelativistic 
constituent ISGW quark model [|1(J. This model is considered to be reliable for a state 
composed of two heavy quarks, which makes it suitable for treating B c ; in addition the 
velocity of B* in the rest frame of B c is to a fair approximation nonrelativistic. In the ISGW 
model the constituent quarks of mass M move under the influence of the effective potential 
V(r) = -4a./(3r) + c + br, c = -0.81 GeV, b = 0.18 GeV 2 0. Instead of the accurate 
solutions of the Schrodinger equation, the variational solutions 



3 3 or —* do. 

ip(r) = 7r~4/?2e ~ or ip{k) = n~^(3~^e 



for S state 



are used, where (3 is employed as the variational parameter. The meson state composed of 
constituent quarks q\ and q 2 is given by 



\M(p))= £ 



1 



2E 



C,sl,s2 



7T 



- - I Mi \M 2 + t 

dkip{k)J — W -^f S 2,si8{p-pi-p2)b[(pi, si, C)d\(p 2 , s 2 , C)|0) 



(13) 

where k is the momentum of the constituents in the meson rest frame, C denotes the colour, 
while f s2yS i = (f I +1 t)/v2 for pseudoscalar and / s2)S i = (f I ~l D/V2, f T, I I for vector 
mesons. Using the normalization of the spinors as in |TJ] , we obtain in the nonrelativistic 
limit 



V(q 2 ) 
Ai(q 2 ) 
Fi(q 2 ) 
Vb c 



F 2 (q 2 



2F 3 (q 2 



U -F 3 (q 2 ) 
2m Bc 



m Bc + m B 
1 + [rn Bc - m B 



1 M b {M c - M u )(3 2 Bc 

M u M c M u m B * u {(3 2 Bc + I3 2 B , 

Fs(q 2 ) , 



m B * ( 2p B J B * 



m B r \0R+0 



J,,, 



( 1 M b (M c + M u ){3 2 Bc 

* J \2M U 2M c M u m B *(f3 2 Bc +f3 2 B * 

2 1 



3M r 3M h 



9b* = m B *— (14) 

and analogously for hb u , /b„ and gs»- Here 



where « = 0.7 flO |. The results for V"(g ) and A\(q ) reproduce the results of [10|, while 
Fi(q 2 ) and F 2 (q 2 ) represent, to our knowledge, the new results within ISGW model. Using 
parameters (3 |L4) and meson masses given in Table 1 and the constituent quark masses 
M u = 0.33 GeV, M c = 1.82 GeV and M b = 5.2 GeV @ we get 

/b„ = 0.18 GeV , flJ = 0.86 GeV 2 , /x B „ = 1.81 GeV' 1 
f Bc = 0.51 GeV , # B * = 2.41 GeV 2 , /x Bc = 0.28 GeV' 1 , 
while the form factors evaluated at q 2 = are given in Table 2. 

2e. The results for 5 C — >■ 5*7 

We use the central value of the current quark masses m u = 0.0035 GeV, m c = 1.25 GeV 
from and V& = 0.04, I4 6 = 0.0035. The SD, VMD and pole contributions to amplitudes 
Ape and Apy needed to compute the amplitude ( ]TD| ) and the decay rate (|TJ) are given in 



Table 3, where the error is due only to the uncertainty in parameter C VMD (|6]). In Ta- 
ble 4 we present the total branching ratio and separately also the SD and LD part of the 
branching ratios for B c — > 5*7 decay, where we have taken t(B c ) = 0.46^°'^ ± 0.03 ps 
as measured by CDF Collaboration recently ||. Note that SD and LD contributions give 
branching ratios of comparable size ~ 10 -8 , which in principle allows to probe the c — ► wy 
transition in B c — > 5*7 decay. Experimental detection of B c — > 5*7 decay at the branching 
ratio well above 10~ 8 would clearly indicate a signal for new physics. The measurement of 
this decay would probe different scenarios of physics beyond the standard model: the non- 
minimal supersymmetric model [|l^] and the standard model with four generations jljj, for 
example, predict 5r(c — > wy) up to 10~ 5 , which would enhance 5r(5 c — > 5*7) up to 10~ 6 . 
In D meson decays (cq — > uq'y), on the other hand, the branching ratios are of order 10 -6 
even within the standard model || ^]: they are driven mainly by the long distance pole 
contributions, which overshadow the c — > wy transition (predicted at the branching ratio 
~ 10 -9 in the standard model) and possible signals of new physics. 

3. B s -> 5^7 DECAY AND s -> d-f TRANSITION 

The calculation of SD and LD contributions to B s — > 5^7 decay is analogous to what 
was presented in the previous section. Due to the small difference mg a — m^*, the final B* d 
meson is almost at rest in the rest frame of B s . 

The SD Lagrangian is 

ZTd 1 = ^ F K^ + m d ^}s F, v , (15) 
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where c^ d7 (m,,) = —0.23 [ Fj . 



The Lagrangian, from which the long distance VMD contribution is calculated, is given 

by 

G 

C-vmd = -^a 2 (m s )V ud V* s d*f(l - 7 5 )s [vry^l - 7s)w - ^(1 - 75 )c] , (16) 
were |a 2 (m s )| ~ 0.5 |l8j and 

g gjg) , 9_M y HM (17) 

VMD 2m2 Qml Y 3m V>* 

From r ea . p (fi~ — > 2^7) < 3.7 • 10~ 9 eV^ the upper limit |CVmd| < 0.01 GeV 2 has been 
obtained JT§. 

The po/e contribution is absent in B s — > £^7 decay since the decay sb — * dfe involves four 
quarks of the same charge. 

We use the constituent quark masses Ma = 0.33 GeV, M s = 0.55 GeV the central 
value of the current quark masses raa = 0.006 GeV, m s = 0.115 GeV [[| and the masses and 
parameters /? of the mesons given in Table 1. Using the formulas of the previous section and 
discarding the pole amplitudes, the resulting amplitudes A PC ,pv and the branching ratios are 
given in Tables 3 and 4, respectively, where the upper limits are due to |CVmd| < 0.01 GeV 2 . 
The upper limit Br(B s — > B^y) < 1.8 • 10~ 7 is dominated by the LD contribution. The same 
conclusion with a smaller upper bound was obtained in |T8j assuming the simple free quark 
decay, which is a reasonable assumption for s —>■ d'y in B s — > B^j decay. However, the cal- 
culation of the VMD amplitude in [|l^] was based on the formalism presented in |19| , which 
is not reliable for s — > dj transition. 



4. SUMMARY 

We have studied flavour changing neutral transitions c — > wy and s — > dj in B c — > B*j 
and B s — > B*^ decays, respectively. The predicted short and long distance contributions to 
these decays within the standard model are presented in Table 4. 

We predict Br(B s — > B^) < 1.8 • 10~ 7 , which is dominated by the long distance contri- 
bution. 

The short distance part (driven by c — > wy) and the long distance part of the branching 
ratio for B c — > B*j decay are found to be of comparable size; they are both of order 10~ 8 . 
This makes B c — > 5*7 decay the most suitable decay to probe the c — > wy transition. Since 
c —>■ wy transition is very sensitive to the physics beyond the standard model, it would be 
very desirable to compare the standard model prediction of Br(B c — > B**y) = (8.5^ 5) • 10~ 9 
presented here to the experimental data in the future. The detection of B c — > B*j decay at 
a branching ratio well above 10~ 8 would signal new physics. In comparison to B c — > B*j 
decay, the D meson decays are far less suitable for probing c — > wy transition, since they are 
almost completely dominated by the long distance effects. 

Finally, we wish to stress that B c — > B*^ and B s — > B*p/ decays are characterized by 
a very clear signature: their detection requires the observation of a B u /Bd decay in coinci- 
dence with two photons. The B c — > 5*7 transition involves the emission of a high energy 
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(985 MeV) and of a low energy (45 MeV) photon in the respective centers of mass of B c , 
B*, while in the decay B s — > £^7 two photons of nearly equal energy (45 MeV) are emitted. 
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Figure caption 

Fig. 1: Long distance contributions in B c — > 5*7 decay, a) VMD contribution; the 
black box denotes the action of the Lagrangian (Bj). b) pole contribution; the white box 
denotes the action of the Lagrangian (§). 
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Table 1: Parameters (3 (taken from ]14|) and masses of pseudoscalar and vector mesons in 
GeV. 
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Table 2: Form factors at q = calculated using ISGW model 10 





ASD 
A PV 


A VMD 
A PV 


a pole 
Apy 


aSD 
A PC 


A VMD 
A PC 


a pole 
A-PC 


B c - B* ul 


5.7 + 17 i 


-14 ± 14 





5.7+17 i 


-7.3 ±7.3 


-21 


B s - B* dl 


62 


< 1.1 • 10 4 





70 


< 5.6 • 10 2 






Table 3: The amplitudes Apv.pc defined in ([TDD f° r SD, VMD and po/e contributions in 
units of 10~ n GeV -1 as predicted by ISGW model. The error-bars in B c — > 5*7 are due to 
the uncertainty in C' VMD = (1.2±1.2) 10~ 3 GeV 2 (|), while the upper bounds for B s — > 5*j7 
are due to |CV.mx>| < 0.01 Ge^ 2 . 
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(8.5+^) -lO" 9 


B s - 


-B* dl 


1.4- KT 11 


< 2.0 • 10~ 7 


< 2.0 • 10" 7 



Table 4: SD, LD and total branching ratios as predicted by ISGW model. The error-bars 
in B c — > fi*7 are due to the uncertainty in C' VMD = (1.2 ± 1.2) 10~ 3 GeV 2 (g), while the 
upper bounds for B s — > B*^ are due to \Cvmd\ < 0.01 GeV 2 . 
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